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ABSTRACT 

A gas- actuated penetration device has been de- 
veloped for high-energy impact testing of structures. 
The high-energy impact testing is for experimen- 
tal simulation of un contained engine failures. The 
non-linear transient finite element code LS-DYNA3D 
has been used in the numerical simulations of a ti- 
tanium rectangular blade with an aluminum target 
plate. Threshold velocities for different combinations 
of pitch and yaw angles of the impactor were ob- 
tained for the impactor-target. test configuration in 
the numerical simulations. Complete penetration of 
the target plate was also simulated numerically. Fi- 
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nally, limited comparison of analytical and experi- 
mental results is presented for complete penetration 
of the target by the impactor. 

INTRODUCTION 

The potential hazard resulting from an un- 
contained turbine engine failure has been a long- 
term concern of the Federal Aviation Administration 
(FA A), National Aeronautical and Space Administra- 
tion (NASA), and the aircraft industry (e.g., Refs. 
[l]-[5]). For the purpose of airplane evaluations, the 
FA A defines an uncontained failure of a turbine en- 
gine as any failure which results in the escape of ro- 
tor fragments from the engine or Auxiliary Power 
Unit (APU) that could result in a hazard (Refs. [4] 
and [5]). A contained failure is one where no frag- 
ments are released through the engine nacelle struc- 
ture; however fragments may be ejected from the en- 
gine air inlet or exhaust. Rotor failures that are of 
concern are those where released fragments have suf- 
ficient energy to create a hazard to the airplane and 
its passengers. 

Accepting that the failures will continue to oc- 
cur in service, attempts are made to contain all debris 
within a strengthened structure (e.g., see Refs. [6]- 
[9] ). Design and test requirements are imposed on the 
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engine nacelle to ensure some containment capabil- 
ity. Engine nacelle design and test requirements are 
covered in the United States Code of Federal Regula- 
tions, Title 14, Aeronautics and Space, Part 33, Air- 
worthiness Standards; Aircraft Engines ([1]). Part 33 
of the Federal Aviation Regulations (FAR) has always 
required the engine nacelle to be designed to contain 
damage resulting from rotor blade failure. The con- 
tainment of failed rotor blades is a complex process 
which involves high energy, high speed interactions of 
numerous locally and remotely located engine compo- 
nents (such as failed blade, other blades, containment 
structure, adjacent cases, bearings, bearing supports, 
shafts, vanes and externally mounted components). 
Once failure begins, secondary events of a random 
nature may occur whose course cannot be precisely 
predicted (e.g., [4]). 

Therefore, assuming that uncontained debris will 
continue to be generated, design considerations out- 
lined in the AC 20-128A (e.g., [5]) provide guidelines 
for achieving the desired objective of minimizing the 
hazard to an airplane from uncontained rotor fail- 
ure. These guidelines assume a rotor failure will oc- 
cur and that analysis of the effects of this failure is 
necessary. The designs intend to make the aircraft 
invulnerable to the debris by such means as deflec- 
tion, the judicious location of critical parts, hydraulic 
lines, and structure, suitable redundency where ap- 
propriate. Given that the damage is uncontained, 
developing an understanding of the impact event of 
the engine fragments or other parts of the structure 
is needed. FAA Advisory Circular AC 20-128A (e.g., 
[5]) provides specifications for fragment sizes to be 
used in the safety analysis models. The fragment size 
includes a single disc with blades fragment, that is, 
with one-third of blade height and one-third of the 
mass of the disc with blades, intermediate fragments 
with one-third of the disc with blades radius with a 
mass of 1 /30th of the disc with blades, and small 
fragments (shrapnel) ranging in size up to a maxi- 
mum dimension corresponding to the tip half of the 
blade airfoil. 

Assuming that the large engine fragments re- 
leased in the radial direction of the engine are con- 
tained by the engine containment structure, exper- 
imental studies and analytical simulations are still 
necessary to understand the effect of small engine 
fragments ejected from the engine air inlet or ex- 
haust on the surrounding structures. Designing air- 
craft structures to either withstand this threat or to 
perform safely after the threat occurs requires an un- 
derstanding of the response of structures subjected 
to high-energy impacts from these small engine frag- 
ments. Although some high-energy penetration work 


has been conducted related to the development of en- 
gine containment structures (e.g., Refs. [l]-[8] ) , there 
is very little reported research (e.g.. Refs. [10]-[12]) 
on metallic and composite airframes when this type 
of high-energy threat occurs. The main objectives 
of the present paper are to discuss the development 
of a high-energy impact device capable of projecting 
small impactor plates that are representative of small 
engine fragments ( 0.6-lb weight), and to present re- 
sults from numerical simulations. The analysis efforts 
presented in this paper address: 

a. The threshold velocity for the impactor to 
penetrate the test specimen when impacted 
at normal incidence to the target. 

b. The threshold velocities for no penetration 
when the impactor strikes the target at at- 
titudes that are combinations of pitch, roll 
and yaw directions. 

c. To compare the analysis results with prelim- 
inary target penetration test results. 

GAS-ACTUATED 
PENETRATION DEVICE 

A photograph of the gas-actuated penetration 
device is shown in Figure 1 (a). The pressure cham- 
ber is connected to the barrel of the device through a 
diaphragm chamber. When the pressure differential 
across the diaphragm reaches a predetermined value, 
the diaphragm ruptures and the gas propels a sabot 
located in the barrel. The sabot is a hollow cylindri- 
cal body made of a plastic material with a provision 
in it to hold an impactor plate. The sabot is guided 
as it travels along the barrel so that the impactor 
plate orientation is not altered significantly before it 
impacts the target test specimen. When the sabot 
reaches the end of the barrel, a splitter arrangement 
in the muffler assembly engages the sabot and re- 
leases the impactor plate which travels farther and 
impacts the target. A photograph of the test speci- 
men mounted in a picture frame fixture in front of the 
barrel is shown in Fig. 1(b). The box surrounding 
the test specimen is used as an impactor containment 
structure and is filled with sand bags. 

The specimen penetration device is designed 
such that the impactor plate located in the sabot is 
projected at the target with specified tolerances on 
its attitude. Since the device will be used for com- 
posite structures which have directional properties, 
these tolerances are necessary to ensure that the im- 
pactor contacts the target, at a given orientation with 
respect to its material axis. The goals for the tol- 
erances in pitch, roll, and yaw angles are 2 degrees 
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which is consistent with the tolerance used in manu- 
facturing laminated composite structures. 

NUMERICAL SIMULATION TOOLS 

Finite element simulations of structural prob- 
lems involve pre-processing, analysis and postpro- 
cessing. For the impact and penetration simulations, 
the pre-processing step is performed using the IN- 
GRID computer code (Ref. [13]). The analysis step 
involves the nonlinear transient dynamic response 
prediction for the dynamic behavior prior to impact, 
the nonlinear impact and penetration event itself, and 
the subsequent dynamic behavior. The LS-DYNA3D 
computer code (Refs. [14, 15]) has been used for the 
analyses here. Additional details are provided in Ref. 
[16]-* 

Penetration of the target plate can be simulated 
in two ways depending on the modeling approach 
used for the target plate. Using the approach of tied 
nodes with failure, coincident nodes are generated in 
selected regions and then tied together with a con- 
straint relation. In the LS-DYNA3D code, these tied 
nodes remain together until the volume-weighted ef- 
fective plastic strain, averaged over all elements con- 
nected to the nodes in a, given constraint, exceeds a 
specified value. This approach is refered to herein as 
the TNWF approach. Once this value is exceeded, 
all nodes in that constraint are released to simulate 
the initiation of a crack, fracture or penetration. In 
the second approach, called the element erosion ap- 
proach, the finite element model is generated in the 
standard manner without requiring duplicate coinci- 
dent nodes. Once the effective plastic strain in an 
element reaches a specified critical value, the element 
is removed from the computations. In this approach, 
elements do not or separate from I he initial finite ele- 
ment model and, hence, tracking the rigid-body mo- 
tion of these newly created fragments is not necessary. 
The LS-TAURUS computer code (e.g.. Refs. [14, 15]) 
is used for postprocessing. 

CONFIGURATION STUDIED 
AND MODELING. 

The target configuration considered in the 
present study is shown in Figure 2. This configu- 
ration is representative of the test setup shown in 
Figure 1 (b). An aluminum target plate is clamped 
between two steel frames and a titanium impactor 
may have pitch and yaw angles as shown in the fig- 
ure. The material properties for aluminum, steel and 
titanium are given in Table 1. An elastic-plastic 
strain-hardening material model is used within LS- 
DYNA3D (Material Type 3 of Refs. [14, 15]). 


The frame that supports the aluminum target 
plate and impactor are modeled using 8-node solid 
elements w r hile the target is modeled using 4-node 
Belytschko-Lin-Tsay shell elements. The node dis- 
tribution in the x, y, z direction is shown in Figure 
2 by numbers in parenthesis. The impactor has a 
node distribution of 9 by 4 by 16 in the x, y, and 
z directions respectively, or 576 nodes with 360 ele- 
ments. Using the INGRID (Ref. [13]) preprocessor, 
the impactor is initially positioned at the center and 
0.25 in. away from the target. The impactor may 
then be pitched or yawned by two degrees as required. 
When the impactor is pitched by two degrees, it has 
to be translated by -0.1362 in. in the v-direction. 
Similarly, when the impactor is yawed by two de- 
grees it has to be translated by -0.1362 in. in the 
z-direction. These translations are necessary since 
the impactor is moving along the its orientation as 
shown in Figure 2. The translations are computed 
based on the distance between the impactor and 
the target in the gas launcher-target, arrangement. 
The velocity of the impactor is Uo cos(a p )cos{(Xy) in 
the x-direct.ion, —Vo sin(a p ) in the y-direction and 
—Vocos(a p )sin(a y ) in the z-direction, where a p and 
a y are the pitch and yaw angles, respectively, and Vo 
is the speed of the impactor. All of the nodes of the 
impactor are prescribed with these velocity compo- 
nents. 

The area with the dashed-line boundary as 
shown in Figure 2 is herein refered to as the shell- 
break area and consists of coincident nodes which 
are tied together with a constraint relation. This ap- 
proach for penetration modeling of tied-nodes-with- 
failure (TNWF) is used to simulate penetration of 
the target by the impactor. In the LS-DYNA3D code 
(Refs. [14, 15]), these tied nodes remain together un- 
til the volume- weighted effective plastic strain, av- 
eraged over all elements connected to the coincident 
nodes in a given constraint, exceeds a specified value. 
The specified plastic strain value for aluminum is 0.2 
which is the ultimate strain at failure. 

According to Figure 2, the element size in the 
shell-break area is 0.05-in. and 0.20-in. in the y- 
direction and z-direction, respectively. The element 
size for the impactor is 0.07267-in. and 0. 33333-in. 
in the y-direction and z-direction, respectively. More 
than four elements in the shell-break area span the 
thickness of the impactor. Therefore, the element size 
in the shell break area is 23% of the smallest dimen- 
sion of the impactor contact surface. A highly refined 
model should have an element size in the contact re- 
gion (or shell-break area) between 20% to 25% of the 
smallest dimension of the impactor contact surface 
based on the studies reported in Ref. [17]. Hence, 
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the model described in Figure 2 is highly refined. The 
finite element model consists of 11,730 nodes, 1,896 
solid elements, and 4,408, shell elements. 

Contact or impact algorithms have always been 
an important capability in the DYNA3D family of 
codes. Contact may occur along surfaces of a sin- 
gle body undergoing large deformation, between two 
or more deformable bodies, or between a deformable 
body and a rigid barrier. In the present study, 
the sliding interface with friction and separation ap- 
proach (LS-DYNA3D, Interface Type 3) is used to 
model the impact event between the impactor and 
target plate, and the friction coefficients are pre- 
scribed to be equal to zero. The bounding surface 
of the three-dimensional impactor is treated as the 
slave surface, and the target plate as the master sur- 
face. 

The LS-DYNA3D code permits automatic ex- 
amination of the finite element mesh and material 
properties in order to determine an appropriate time 
step size for numerical stability. This time step size is 
then automatically adjusted throughout the transient 
analysis to account for contact and local material and 
geometric nonlinearities. 

RESULTS AND DISCUSSION 

Analytical simulation results obtained using the 
LS-DYNA3D code are reported in this section for the 
target configuration described in Figure 2. The mod- 
eling features discussed in the previous section were 
used. The finite element studies are based on the as- 
sessment of the time variation of the axial velocity 
of the centroid of the impactor, the time variation 
of the contact force in the axial direction, and the 
maximum plastic strain on the target. 

Results are presented for the following cases with 
different initial impactor velocity (Vo) and different- 
pitch (Op) and yaw (ay) angles of the impactor; 

• Case 1: Vo = 5,400 in/sec (450 ft /sec), a p = 0 
degrees, a y = 0 degrees. 

• Case 2: Vo = 3,000 in/sec (250 ft/sec), a p — 0 
degrees, a y ~ 0 degrees. 

• Case 3: Vo = 3,000 in/sec (250 ft/sec), ot p = 0 
degrees, a y ~ 2 degrees. 

• Case 4: Vo = 2,640 in/sec (220 ft/sec), a p = 2 
degrees, a y ~ 0 degrees. 

• Case 5: V'o = 2,760 in/sec (230 ft/sec), a p = 2 
degrees, a y = 2 degrees. 

Case 1 is intended to simulate the penetration 
test with the aluminum target plate. The damage 
result from this analysis case are used to compare 
with the corresponding experimental results. 

Analytical simulations for Case 2 through 5 were 
performed to determine the threshold velocity for 


penetration for the different pitch and yaw angle com- 
binations. The threshold velocity is defined as the 
velocity above which the impactor will penetrate the 
target completely. The threshold velocity is an im- 
portant quantity in selecting the initial impactor ve- 
locity for the gas launcher. The threshold velocity 
was determined by analyses to be the velocity for 
which the impactor rebounds from the target while 
creating partial perforation of the target or produc- 
ing a maximum plastic strain that is marginally lower 
than the ultimate strain value of 0.2 for aluminum 
without penetration of the target. 

Simulation Parameters 

To insure that a simulation analyses using the 
LS-DYNA3D code is meaningful, the analyst needs 
to monitor the time step size, the ratio of the sliding 
interface energy to the initial or total energy, and the 
nodal velocities. Since an explicit time integration 
algorithm is used in this code which automatically 
adapts the time step size as the plasticity and dam- 
age develop in the elements, the time step size may be 
driven to nearly zero. A problem with the simulation 
occurs for such a case. Similarly, if the ratio of the 
sliding interface energy to the total energy (Max(SIE 
/ TE)) is larger than 10%, then the sliding interface 
penalty factor (SIPF) has to be adjusted. A good 
goal is to keep this ratio to be under 10%. A high ra- 
tio of the sliding interface energy to the total energy 
may lead to a simulation problem where nodal veloc- 
ities take on out-of-range values (e.g., Not a Number 
values). To achieve these goals in simulation, mul- 
tiple analyses with the current finite element model 
and possible finite element remodeling is necessary to 
validate the simulation results. 

The sliding interface penalty factor (SIPF) and 
the time step scale factor (TSSF) are two parameters 
among others that affect the simulation results. The 
variation in simulation results for Case 5 for different 
SIPF values is shown in Table 2. The value for TSSF 
is 0.6 for all analyses. This value for TSSF does not 
lead to spikes in the time variation of the impactor 
axial velocity. Hence TSSF — 0.6 is an appropriate 
value for the impactor velocities considered in Table 
2. This value is also appropriate for the impactor 
velocities considered in the other simulations. 

It can be seen from Table 2, that there are no 
simulation problems with Analysis Number 1 and 2, 
but the Max(SIE / TE) values are 33% and 12% re- 
spectively. For Analysis Number 3 and 4, although 
the Max(SIE / TE) values are well below 10%, there 
are simulation problems. Analysis Number 5 and 6 
were carried out for a velocity of 230 ft./sec and al- 
though the Max(SIE / TE) value is 4% for Analy- 
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sis Number 5, there is a simulation problem. The 
Max(SIE / TE) value is 8.5% for Analysis Number 
6 which is considered to be acceptable and there is 
no simulation problem. In Analysis Number 6, the 
impactor rebounds from the target while creating a 
maximum plastic strain of 0.1972. Hence Analysis 
Number 6 is considered to provide meaningful results 
for the case when a p = a y = 2 degrees. Since the 
maximum plastic strain is 0.1972 for Analysis Num- 
ber 6, Vo = 230 ft /sec is considered close to the 
threshold velocity. 

Simulation for other cases were investigated in a 
similar manner t o Case 5. The simulation parameters 
and results for Cases 1 through 5 are summarized in 
Table 3. Reference [17] shows the effect of SIPF on 
the time variation of the total energy. Accordingly, 
there is a discontinuity in the total energy associated 
with the out-of-range velocities. The importance of 
simulation parameters for velocities that are close to 
the threshold velocity is also discussed in Ref. [17]. 
It can be seen from Table 3 that the Max(SIE / TE) 
value is less than 10% for the cases considered. No 
simulation problems occurred up to the analysis ter- 
mination time. At the analysis termination time, the 
impactor is away from the target and there is no con- 
tact with the targft at all. 

Simulation Results 

For Case 1 . when the velocity of the impactor is 
450 ft/sec and with zero pitch or roll, there is com- 
plete penetration of the target and the residual ve- 
locity of the impact or is 334.7 ft./sec. The impactor 
rebounds from tin* target for Cases 2 through 5 with 
residual velocities of -13.6, -99.5, -77.0, and -89.8 
ft/sec, respectively. For Case 2, there is partial pen- 
etration of the target, whereas for Cases 3, 4 and 5, 
there is no partial penetration. However, in Cases 
4 and 5, the maximum effective plastic strains are 
0.1950 and 0.ID72. respectively, which is close to the 
ultimate strain value of 0.2. Therefore, the initial ve- 
locities for Case 1 and 5 are considered close to the 
threshold velocity |or the pitch and yaw angles con- 
sidered. The initial \e|ocjf.y for Case 3 is not close to 
the threshold \e|oeitv for the yaw angle considered, 
since the maximum effective plastic strain is well be- 
low the ultimat* ni ram of 0.2. 

The axial \e|ocii\ histories for the centroid of the 
impactor for ( *as»^ I 1 1 trough 5 are shown in Figure 3 
and the axial o>nia< t force for the target for Cases 1 
through 5 are shown in Figure 4. The residual veloc- 
ity for Case 1 is positive and indicates penetration, 
while for the other rases the residual velocities are 
negative indicating rebound from the target. The dif- 
ference in the velocity history between Case 1 and the 


other cases is mainly due to the large difference in the 
velocity of the impactor which leads to penetration of 
the target. The differences between the axial veloc- 
ity histories for Cases 2 through 5 are small due to 
the small differences in impactor velocity and pitch 
and yaw angles. The maximum value of the axial 
contact force for Case 1 is much greater than for the 
other cases. The duration of contact is much shorter 
than for other cases, since the impactor velocity for 
Case 1 is greater than for the other cases. The small 
differences in axial contact force histories for Cases 2 
through 5 are due to differences in impactor velocities 
and pitch and yaw angles. 

A close-up view of the target deformed geome- 
try (at t-— 1 ms) in the vicinity of the impact site for 
Case 1 is shown in Figure 5, where petaling at the 
backside of the target can be seen. A contour plot 
of the effective plastic strain results (at t— 1 ms) for 
Case 1 is presented in Figure 6. A contour plot of the 
effective plastic strain results (t=1.8 ms) for Case 2 
is presented in Figure 7, which shows the partial pen- 
etration of the target. 

A test was conducted with an impactor axial ve- 
locity of 450 ft/sec, and the impactor penetrated the 
aluminum target. On examining the damaged target, 
it was observed that the impactor made contact with 
the target below the center of the target. This ob- 
servation suggests that the impactor had a different 
pitch angle. The exact pitch angle is difficult to assess 
since some test conditions were unknown. A photo- 
graph of the region around the impact site of the test 
target is shown in Figure 8. When comparing Figure 
8 with Figure 5, it can be seen that the deformed ge- 
ometry for Case 1 (Figure 5), which has zero pitch 
and yaw angles, has less petaling on one side of the 
impact site in Figure 5 than in Figure 8. The dif- 
ference in damage shown in Figure 8 may be due to 
the impactor having a pitch angle which is not small. 
Also comparing the simulation results for Cases 2 and 
4, it can be seen that a two-degree pitch angle in Case 
4 led to a decrease from 250 ft/sec (threshold velocity 
for Case 2 ) to 230 ft /sec (threshold velocity for Case 
4 ) or a 12% decrease in threshold velocity between 
Cases 2 and 4. This obsevation suggests that the test 
configuration is sensitive to even small pitch angles 
for the impactor. 

CONCLUDING REMARKS 

A gas- actuated penetration device has been de- 
veloped for high-energy impact testing of structures. 
High-energy impact tests were conducted to deter- 
mine the threshold velocities for complete penetra- 
tion of the target plate by a rectangular titanium 
plate impactor. Threshold velocities for different 
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combinations of pitch and yaw angles of the impactor 
were obtained for the impactor-target test configu- 
ration in the numerical simulations. The numerical 
simulation results indicate that the threshold velocity 
is sensitive to small pitch and yaw angles of the im- 
pactor. A pitch angle of the impactor decreases the 
threshold velocity compared to the threshold veloc- 
ity for zero pitch and zero yaw angle of the impactor. 
A yaw angle of the impactor increases the threshold 
velocity compared to the threshold velocity for zero 
pitch and zero yaw angle of the impactor. Numerical 
simulation predicted damage similar to that obtained 
from an experiment for complete penetration of the 
target by the impactor. 
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Table 1 Material properties for aluminum, titanium, and steel. 



Aluminum 

(2024-T3) 

Titanium 

(T1-6A1-4V) 

Steel 

(A36) 

Young’s modulus (Msi) 

10.0 

16.0 

30.0 

Poisson’s ratio 

0.30 

0.30 

0.30 

Yield strength (Msi) 

0.05 

0.12 

0.35 

Tangent modulus (Msi) 

0.10 

0.30 

0.15 

Hardening parameter 

0.20 

0.20 

0.23 

Weight density (lbs/in 3 ) 

0.10 

0.16 

0.28 

Ultimate strain to failure (in/ in) 

0.20 

- 

- 


Analysis 

No. 

1 

2 


3 


4 

5 

6 


Table 2 Effect of SIPF on numerical simulation results for Case 5. 

Initial velocity, SIPF a TSSF* Max(SIE / TE) C Termination 

Vo, ft /sec time, ms 

250 0.25E-03 0.60 33% 1.8 

Remark: Maximum effective plastic strain is 0.152, impactor rebound. 

250 0.35E-03 0.60 12% 1.8 

Remark: Maximum effective plastic strain is 0.187, impactor rebound. 

250 0.55E-03 0.60 4.8% 1.8 

Remark: Target plate penetration and out-of-range velocities at 1.5 ms. 

250 0.45E-03 0.60 6.1% 1.8 

Remark: Target plate penetration and time step size driven to zero at t— 1.317 ms. 

230 0.45E-03 0.60 4%. 1.8 

Remark: Target plate perforation and impactor rebound. Out-of-range velocities at 1.14 ms 
230 0.35E-03 0.60 8.5% 1.9 

Remark: Maximum effective plastic strain is 0.1972 and impactor rebound. 

a SIPF = sliding interface penalty factor 
& TSSF = time step scale factor 

c Max(SIE / TE ) = maximum sliding interface energy to total energy 


Table 3 Summary of simulation parameters and results for Cases 1 through 5. 


Case No. 

1 

2 

3 

4 

5 

Initial velocity, Vo (ft/sec) 

450 

250 

250 

220 

230 

(a p ,a y ) degree 

(0,0) 

(0,0) 

(0,2) 

(2,0) 

(2,2) 

Termination time (ms) 

1.0 

1.8 

1.8 

1.8 

1.9 

SIPF 

0.5E-03 

0.75E-03 

0.25E-03 

0.65E-03 

0.35E-03 

TSSF 

0.6 

0.6 

0.6 

0.6 

0.6 

Max(SIE / TE) % 

1.0 

2.2 

2.2 

6.0 

8.5 

Residual velocity (ft/sec) 

334.7 

-43.6 

-99.5 

-77.0 

-89.8 

(axial direction) 
Maximum plastic strain 

0.2000 

0.2000 

0.1460 

0.1950 

0.1972 

Target penetration 

Complete 

partial 

NO 

NO 

NO 

Impactor rebound 

NO 

YES 

YES 

YES 

YES 

Number of cycles 

8308 

14954 

14954 

14954 

15785 

CPU* time (sec) 

♦ riTITT 

4591 

2 . ’ • r nr 

8104 

'IT T 1* TV 

8272 

8210 

8421 


*CPU time is for SGI Indigo R4400 machine 
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(a) Projectile launcher 
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(b) Target plate assembly 


Figure 1: Photograph of the gas-actuated projectile launcher and test setup. 
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Figure 2: Dimensions and finite element node distribution for the test configuration. (Not to scale. All 

dimensions in inches.) 
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Time, ms 


Figure 3: Axial velocity history of impactor for simulation Cases 1 through 5. 




Figure 5: Deformed geometry of target in the vicinity of the impact site for simulation Case 1 at t = 
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Figure 6: Contour plot of effective plastic strain for simulation Case 1 at t = 1 ms 
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Figure cS; Photograph of damage aluminum target for an impactor initial velocity of 450 ft/s 
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